Analysis of multi-phase mixtures

ABSTRACT

A method for continuous analysis of mixtures for determining selective properties of components in the mixture flowing along a flow path includes the steps of, firstly, directly measuring and/or calculating the pressure, temperature, density and mass flow rate of the mixture at a first point upstream of a characteristic change zone in the flow path. A change in the characteristics of the mixture flow is caused in the characteristic change zone, and thereafter, the pressure, temperature, density and mass flow rate are measured directly and/or calculated at a second point downstream of a characteristic change zone. Selected properties of the mixture are determined on the basis of the information collected together with known parameters using algorithms.

BACKGROUND AND FIELD OF THE INVENTION

This invention relates generally to a method and apparatus for determining selected properties of a multi-phase mixture.

Whilst it will be appreciated by those persons skilled in the art that the method and apparatus of the invention are readily applicable to the analysis of a wide variety of multi-phase mixtures, one particularly advantageous application relates to the analysis of the oil mixture discharged from an oil well. It will be convenient therefore to hereinafter describe the invention with reference to that particular application.

It is often necessary to be able to monitor and/or determine the various properties of the oil/gas/water mixture being discharged from an oil well. The properties of the mixture which are particularly important are the volume rate and density of each phase. It is the current practice when analysing the mixture discharged from an oil well to feed a relatively large quantity of the mixture to a "test separator". The "test separator" separates the multi-phase mixture into discrete phases before flow measurement is effected.

Separators of this general type are known and it is not proposed to provide a detailed description here. The major disadvantages of separators of this type are that they are relatively bulky and take up considerable space particularly when used on off-shore platforms. As a result each offshore platform usually has only one test separator, and this test separator is employed to take measurements from one well at a time. A further disadvantage is that they are relatively expensive.

It is an object of the present invention to provide a relatively simple method and apparatus for analysing multi-phase mixtures and with reference to the specific application referred to above, it is an object to alleviate one or more of the aforementioned disadvantages.

SUMMARY OF THE INVENTION

According to one aspect of the present invention there is provided a method of analysing multi-phase mixtures for determining selected properties of components in the mixture flowing along a flow path comprising the steps of directly measuring or calculating the density of the mixture at two or more points along the mixture flow path, directly measuring or calculating the mass flow rate of the mixture at at least one of the aforementioned points along said flow path and causing a change in the characteristics of the mixture flow between said one or more adjacent points and thereafter determining the selected properties on the basis of information collected.

The characteristic which is changed can be any one of a suitable group. For example the characteristic changed may be a pressure change or a temperature change or at least a partial separation of the phases in the mixture.

The characteristic change may also be effected by adding a substance to the mixture between the two points. This substance may be in the form of a chemical to cause a reaction or may be by the provision of a side stream of gas. In another arrangement, the characteristic change can be initiated by a change at some other point in the flow path which thereby causes the characteristic change. For example, a pressure differential upstream of the two points may be sufficient to cause a change between the points.

By using the above method, it is possible to determine the properties of the various phases of the mixture such as, for example, the volume and density of each particular phase. To determine these properties it may be necessary to carry out more than two density and/or mass flow measurements or calculations and where two or more measurements are carried out a corresponding flow characteristic changes are required between the adjacent measuring points. In some instances not all variables may be required to be determined purely from the mass flow measurement technique. For example, in mixtures such as the oil mixture from oil wells it may be convenient to assume that the density of the water phase and the oil phase remain constant along each pressure differential. From the information collected, it may be possible to produce suitable algorithms from which the required characteristics can be determined.

Thus in the particular example of a mixture discharged from an oil well, a series of equations can be written from which the various properties of the phases of the mixture can be determined.

In a typical example, it is reasonable to assume that the density of the oil and water remain substantially constant throughout the mixture flow. As such, it is possible to determine these values in any suitable manner. Using the method of the present invention the following properties are measured:

(a) The density of the mixture at points 1 and 2 (ρ₁ and ρ₂).

(b) The mass flow rate of the mixture (M_(T1)) (where the characteristic change is a pressure differential, the mass flow rate remains constant).

(c) The temperature and pressure at points 1 and 2.

(d) The density of the oil (ρ_(o) and ) and water (ρ_(w)) (As mentioned earlier the density of the oil and water are measured separately).

    P.sub.1 =k.sub.1 ρ.sub.1 T.sub.1                       ( 1)

    P.sub.2 =k.sub.2 ρ.sub.2 T.sub.2                       ( 2)

    M.sub.T1 =M.sub.T2                                         ( 3)

    M.sub.T1 =ρ.sub.g1 V.sub.g1 +ρ.sub.w1 V.sub.wl +ρ.sub.o1 V.sub.ol                                                  ( 4)

    M.sub.T2 =ρ.sub.g2 V.sub.g2 +ρ.sub.w2 V.sub.w2 +ρ.sub.o2 V.sub.o2                                                  ( 5)

    V.sub.w1 =V.sub.w2                                         ( 6)

    V.sub.o1 =V.sub.o2                                         ( 7)

    V.sub.T1 =V.sub.g1 +V.sub.o1 +V.sub.w                      ( 8)

    V.sub.T2 =V.sub.g2 +V.sub.o2 +V.sub.w2                     ( 9)

V is volume

subscript:

g is gas

o is oil

w is water

If ρ_(w1), ρ_(o1), ρ_(w2), ρ_(o2), k₁, and k₂ are assumed or hand measured it is therefore possible to determine V_(g1), V_(g2), ρ_(g1), ρ_(g2), V_(o1), V_(o2), V_(w1), V_(w2), V_(T1), V_(t2).

According to another aspect of the present invention there is provided apparatus for analysing multi-phase mixtures comprising means for directly measuring the density of a mixture at at least two selected points along a flow path of the mixture, means for measuring the mass flow rate of the mixture at at least one of the points and means for using a change in the flow characteristics of the mixture between the points.

The apparatus of the present invention may comprise means for directly measuring the density of a mixture at at least two selected points along the path of travel of the mixture, means for measuring the mass flow rate of the mixture at at least one of the points and means for causing a change in the flow characteristics of the mixture between the points. It will be appreciated that this apparatus can be interfaced with a suitably programmed computer so that, for example, the properties of the mixture utilizing the equations above can be quickly determined.

The means for directly measuring the density and mass flow rate of the mixture may be in the form of a mass flow meter. A typical example of a suitable mass flow meter is sold under the trademark Micro Motion and is described in U.S. patent specification No. 422,338, now abandoned, and U.S. Pat. No. Re. 31,450. It may also be necessary to measure the temperature and pressure of the mixture at the mass flow meters and this can be done in any suitable manner. In another arrangement the density of the mixture could be measured by an on-line density meter one example of which is sold under the trademark SOLARTRON by the company Solartron Schlumberger.

If the characteristic change is a pressure differential then this can be caused by any suitable device such as a choke or the like.

If the characteristic change is a phase separation then in one preferred form a hydro cyclone of the type which is for example described in International application PCT/AU84/00195 which has been specifically developed for the treatment of oily water can be used. In one form the mass flow measurement may be taken at the inlet to the hydro cyclone and the outlet therefrom containing the oil/gas phase which is removed from the oil water mixture. In certain instances it may be necessary to use additional mass flow measurement meters with two or more hydro cyclones in series.

In the above described arrangements the mass flow measurement is of the total flow. It is possible however to merely take a sample of the flow and pass that sample through a device comprising two mass flow measurement meters with a pressure differential therebetween. Using such an arrangement however it is only possible to determine the relative percentage volumes of the phases.

BRIEF DESCRIPTION OF THE DRAWINGS

Preferred embodiments of the invention will hereinafter be described with reference to the accompanying drawings in which:

FIG. 1 is a schematic circuit drawing showing the general principle of the invention;

FIG. 2 is a schematic circuit drawing incorporating one form of device according to the invention;

FIG. 3 shows an arrangement by which a sample of the mixture is taken from the main line;

FIGS. 4 and 5 show a form of the apparatus when used with hydro cyclones;

FIG. 6 is a circuit diagram showing a use of the apparatus in relation to gas optimization;

FIG. 7 is a circuit diagram showing a use of apparatus for controlling discharge from an oil well head; and

FIG. 8 is a circuit diagram of apparatus used in the example.

DETAILED DESCRIPTION OF THE PREFERRED EMOBODIMENTS

FIG. 1 shows schematically what may be required to be measured on-line at points 1 and 2 and with a characteristic change and C to determine the desired properties of the various phases in the mixture.

Referring to FIG. 2 of the drawings, there is shown an oil rig 1 comprising oil wells 2, 3 and 4, a discharge line 6 for delivering the mixture to a separator (not shown). The device of the present invention comprises a pair of mass flow meters 10 and 11 having a choke 13 therebetween. The mass flow meters can measure the mass flow as well as the density of the mixture.

In FIG. 3, there is shown an embodiment where a sample is taken from the main delivery line 6, the sample passing through mass flow meters 10 and 11 and choke 13.

In FIGS. 4 and 5 there is shown the arrangement of mass flow meters 10, 11 13 and 14 with hydro cyclones 15 and 16. In this particular arrangement the characteristic change as effected by phase separation. For example, in FIG. 4 mass flow meter 10 measures the mass flow and density of the mixture in total and mass flow meter 11 measures the mass flow and density of the oil and gas phases. It is therefore possible to device equations to determine the desired properties of each phase.

One application of the invention relates to gas optimization of oil wells. To bring the oil water mixture to the surface of a well, gas is delivered down the well and thereby reduces the mass of the mass of the column of mixture within the production tubing. As a result of this mass reduction, the reservoir pressure causes the column to rise to the surface.

The amount of gas fed down the well will affect the oil/water/gas mixture delivered to the surface. Thus, by controlling the gas flow the desired oil/water/gas mixture can be controlled. FIG. 6 is one example of how the gas flow can be controlled using the method and apparatus of the presnet invention.

Gas is fed under pressure along line 30 through throttling valve 20 down the well where it mixes with the column of oil/water and sometimes gas to bring the mixture to the surface. The properties of the mixture are determined through mass flow meters 10 and 11 and choke valve 13 and the information fed to controller C. The controller then opens or closes throttling valve 20 as required.

In FIG. 7 there is shown an arrangement whereby well flow is controlled by the control of the well head choke. In this installation, two mass flow meters 10 and 11 are operatively connected to the well output line, the meters being disposed on opposite sides of the choke 13. The instrumentation provided at each point consists of a temperature transmitter, a pressure transmitter and a density transmitter.

The transmitter measurements are adapted to be fed to a processor P such as a computer containing data and algorithms regarding characterization of gas pressure/temperature/volume/density relationships, mass balance relationships and data regarding oil and water properties as a function of temperature and pressure. The processor P uses this data and the transmitter signals to determine oil, water and gas mass and volumetric flow rates. This information is used as a basis for controlling the position of the choke. This may be done automatically or manually.

By employing the multi-phase analysis system of the invention, for example, on an off-shore platform with multiple wells, continuous measurement of production from each well can be achieved. This can lead to improved reservoir management. The analysis also provides for a more representative flow measurement from the well head due to an improved response time compared to the "test separator" discussed earlier.

In the particular example given above, the system is likely to be able to monitor slugging flow whereas the "test separator" cannot effectively owing to its poor response time.

An example of multi-phase measurement according to the present invention is given below for measurement of a distillate fuel oil/water/air mixture, the measurement being done in a laboratory.

FIG. 8 shows the layout of the test installation. The mixture tested was a mixture of diesel, water and air. The water, diesel and air were fed into the system via flow meters M₁, M₂ and M₃ and through gate valves V₈, V₉ and V₁₀ respectively. M₁ was a Fisher and Porter Rotameter M₂ a VAF positive displacement vane meter and M₃ a Fisher and Porter Purge Rotameter. The gate valves were used to vary the inlet flows.

The flow of the various phases passed to a static flow mixer 1 which provided for a homogeneous three phase flow. A clear pipe section 2 was provided for enabling visual inspection.

The mixture passed to the test circuit section of the installation where measurement of the flow rate was made by meter M₄ which was a Micrometer D40 mass flow meter with the DT7 density and temperature and DT10 flow rate and flow stabilization units. The flow meters M₁, M₂ and M₃ were calibrated against meter M₄ before the test and the meter M₄ was volumetrically calibrated before the test.

The pressure gauge ρg measured the pressure at the meter M₄ and the pressure gauge ρa was used in conjunction with meter M₃ to give the inlet air flow.

The valves V₁, V₂, V₃, V₄, V₅ and V₆ were arranged to control the direction of flow and V₇ was the pressure dropping gate valve. This valving arrangement was used to simulate upstream and downstream of the pressure dropping valve V₇. As such, it was only necessary to use one meter M₄ rather than two.

Initial readings were taken of the three inlet flows and upstream conditions at the meter M₄ thereafter the system was reversed to simulate downstream conditions and further readings were taken. Table 1 sets out the results of four tests conducted.

The mathematics involved with test number 2 are given below.

Initial Known Data

Density of water (ρw): 1.00s.g.

Density of diesel (ρd): 0.828s.g.

Density of air (ρg): 1.29×10⁻³ (0° C. 1 Atmos)

    ______________________________________                                                              After Valve V.sub.7                                       Before Valve V.sub.7 (UPSTREAM)                                                                     (DOWNSTREAM)                                              ______________________________________                                         M.sub.T1                                                                             19.2        Kg/min     M.sub.T2                                                                              19.2                                       P.sub.1                                                                              330         KPng       P.sub.2                                                                               265.                                       T.sub.1                                                                              19.4        °C. T.sub.2                                                                               19.4                                       ρ.sub.T1                                                                         .825        s.g.       ρ.sub.T2                                                                          .812                                       ______________________________________                                    

Inlet Conditions

V_(d) '10.7 l/min:V_(w) ' 10.4 l/min:V_(g) ' 2.2 l/min

The above flows were corrected to calibrated flows from original readings.

Calculating volume flows from above ##EQU1## ΔV_(T) =V_(T2) -V_(T1) =0.4 l/min V_(T1) =V_(d) '+V_(w) +V_(g) '=23.3 l/min

The expected change in V_(T) due to the change in gas volume induced across V_(T) can be calculated from PV=n RT (ideal gas).

For this case n RT is constant. ##EQU2##

(This reuult compares with ΔV_(T) as measured by M₄ of 0.4 l/min).

Given that the volume flow rates of diesel and water will remain constant across V₇ ; that is

    V.sub.d1 +V.sub.w1 =V.sub.d2 +V.sub.w2

then the measured change in volume flow rate across V₇ is the flow rate change in the gas phase; that is ##EQU3## Using the ideal gas law in this instance ##EQU4## From the above result ΔV_(g) =0.4 l/min

T1=T2=19.4° C.

P1=330 KPag

P2=265 KPag

nR=3.3

At meter M₃ Pa=470 KPag T=20° C. ##EQU5## This compares with the calibrated inlet volume of 1.7 l/min. Changing the above flow rate to the flow rates at the upstream and downstream

V_(g1) =2.2 l/min

V_(g2) =2.6 l/min

To find the gas density at the above conditions it is necessary to know its molecular weight

    M=29 kg/kmol

    ρ.sub.g1 =5.1 g/l

    ρ.sub.g2 =4.3 g/l

The calculation of the water and diesel flow rates can be found from the following equations ##EQU6## Substituting values into the above

    V.sub.w =10.5

    V.sub.d =10.5

These are comparable to the inlet flows.

                                      TABLE 1                                      __________________________________________________________________________                         AFTER                                                      BEFORE VALVE (UPSTREAM)                                                                            VALVE (DOWNSTREAM)                                                                              INLET CONDITIONS                            M.sub.T1                                                                           P.sub.1                                                                           T.sub.1                                                                           Density         Density  Diesel                                                                             Water                                                                              Air                                 kg/min                                                                             kPag                                                                              °C.                                                                        ρ.sub.T1                                                                        V.sub.T1                                                                          M.sub.T2                                                                          P.sub.2                                                                           T.sub.2                                                                          ρ.sub.T2                                                                         V.sub.T2                                                                          V.sub.d '                                                                          V.sub.w '                                                                          V.sub.g '                                                                         V.sub.T1 '                                                                        M.sub.T1 ˜                                                               ρ.sub.d V.sub.d '                                                          + ρ.sub.w V.sub.w       __________________________________________________________________________                                                        '                           (1)                                                                              19.0                                                                               325                                                                               19.5                                                                              .873 21.8                                                                              19.0                                                                              260                                                                               19.4                                                                             .863  22.0                                                                              10.8                                                                               9.9 .91                                                                               21.6                                                                              ˜18.7                 (2)                                                                              19.2                                                                               330   .825 23.2                                                                              19.2                                                                              265  .812  23.6                                                                              10.7                                                                               10.4                                                                               1.72                                                                              22.8                                                                              ˜19.3                 (3)                                                                              19.0                                                                               360   .775 24.3                                                                              19.0                                                                              295  .759  25.1                                                                              10.8                                                                               9.8 3.7                                                                               24.4                                                                              ˜18.7                 (4)                                                                              17.4                                                                               325   .751 23.2                                                                              17.4                                                                              215  .719  24.3                                                                              14.5                                                                               5.2 3.4                                                                               23.1                                                                              ˜17.2                 __________________________________________________________________________                           ##STR1##                                                                              ΔV.sub.T = V.sub.T2 - V.sub.T1                                                     and pressures at measurement                                                  points.Δvg = V'gl - Vg2,                                                 Vgl,2 calculated from gas                __________________________________________________________________________                                           law                                                         (1)      0.2       0.2                                                         (2)      0.4       0.4                                                         (3)      0.6       0.6                                                         (4)      1.1       1.2                                      __________________________________________________________________________ 

We claim:
 1. A method for continuous analysis of mixtures for determining selective properties of components in the mixture flowing along a flow path, the method including the steps of:(a) Firstly, directly measuring and/or calculating the pressure, temperature, density and mass flow rate of the mixture at a first point upstream of a characteristic change zone in the flow path; (b) causing a change in the characterisitcs of the mixture flow in said characteristic change zone; (c) thereafter measuring directly and/or calculating the pressure, temperature, density and mass flow rate at a second point downstream of said characteristic change zone; and (d) determining the selective properties of the mixture on the basis of the information collected together with known parameters using algorithms.
 2. A method according to claim 1 wherein the characteristic which is changed is a pressure change.
 3. A method according to claim 1 wherein the characteristic which is changed is a temperature change.
 4. A method according to claim 1 wherein the characteristic which is changed is at least a partial separation of the phases in the mixture.
 5. A method according to claim 1 wherein said characteristic change is initiated by a change other than between said two points which causes said characteristic change.
 6. A method according to claim 1 including the step of utilizing an interfaced computer to determine the selected properties.
 7. Apparatus for continuous analysis of mixtures for determining selective properties of components in the mixture flowing along a flow path, the apparatus comprising:a first mass flow meter for directly measuring the density and mass flow rate at a first point upstream of a characteristic change zone in the flow path; means for causing a change in the characterisitcs of the mixture flowing in said characteristic change zone; and a second mass flow meter for directly measuring the density and the mass flow rate of the mixture at a second point downstream of the characteristic change zone, the arrangement being such that selected properties of the mixture can be determined on the basis of the information collected together with known parameters using algorithms.
 8. Apparatus according to claim 7 including a computer arranged so that the properties of the mixture can be determined from the information obtained from the other parts the apparatus.
 9. Apparatus according to claim 8 wherein said means for directly measuring the density and mass flow rate of the mixture is in the form of a mass flow meter (10, 11).
 10. Apparatus according to claim 7 wherein said means for directly measuring the density and mass flow rate of the mixture is in the form of a mass flow meter (10, 11).
 11. Apparatus according to claim 7, wherein the characteristic changed is a pressure differential caused by a choke valve (13).
 12. Apparatus according to claim 7 wherein the characteristic changed is a phase separation by a hydro cyclone (15, 16).
 13. A method of controlling discharge from a well head comprising the steps of analysing the multi-phase mixture being discharged from the well head in accordance with claim 1 and from the results obtained controlling the discharge conditions in a selected manner.
 14. A method according to claim 13 wherein the discharge of the mixture from the well is controlled by a choke, the method including operating the choke in response to the results obtained.
 15. A method according to claim 13 wherein the discharge of the mixture from the well is controlled by the delivery of gas to the well the method including controlling the delivery of the gas in response to the results obtained.
 16. Apparatus for controlling discharge from a well head the apparatus being characterized by multi-phase analysis apparatus for analysing multi-phase mixtures according to claim 7 a choke (13) operatively connected to the discharge line from the well head and control means (P) for controlling opening and closing of said choke said control means being operable in response to information received from said multi-phase analysis apparatus.
 17. Apparatus for controlling discharge from a well head, the apparatus being characterized by multi-phase analysis apparatus according to claim 7 means (30) for delivering gas to the well having valve means (20) controlling the quantity of gas delivered and control means (C) for controlling opening or closing of said valve means said control means being operable in response to information received from said multi-phase analysis apparatus. 